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Baryon mass extrapolation 
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Consideration of the analytical properties of pion-induced baryon self-energies leads to new functional forms 
for the extrapolation of light baryon masses. These functional forms reproduce the leading non-analytic behavior 
of chiral perturbation theory, the correct heavy-quark limit and have the advantage of containing information on 
the extended structure of hadrons. The forms involve only three unknown parameters which may be optimized by 
fitting to present lattice data. Recent dynamical fermion results from CP-PACS and UK-QCD are extrapolated 
using these new functional forms. We also use these functions to probe the limit of the chiral perturbative regime 
and shed light on the applicability of chiral perturbation theory to the extrapolation of present lattice QCD 
results. 
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1. FORMALISM 

In recent years there has been tremendous 
progress in the computation of baryon masses 
within lattice QCD. Still, it remains necessary to 
extrapolate the calculated results to the physi- 
cal pion mass (fj, = 140 MeV) in order to make 
a comparison with experimental data. In doing 
so one necessarily encounters some non-linearity 
in the quark mass (or m^), including the non- 
analytic behavior associated with dynamical chi- 
ral symmetry breaking. We recently investigated 
this problem for the case of the nucleon magnetic 
moments [Q. It is vital to develop a sound un- 
derstanding of how to extrapolate to the physical 
pion mass. 

1.1. Self-Energy Contributions 

Chiral symmetry is dynamically broken in 
QCD and the pion alone is a near Goldstone bo- 
son. It is strongly coupled to baryons and plays 
a significant role in N and A self-energies. The 
one-loop pion induced self-energies of the N and 
A are given by the processes shown in Fig. |l|. 

In the standard heavy baryon limit, the ana- 
lytic expression for the pion cloud correction to 
the masses of the N and A are of the form iHl 
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CiVAf + CTJVA ; M/a = CTAA + CTAJV , (1) 
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Figure 1. One- loop pion induced self-energy of 
the nucleon and the delta. 
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Here AM = Ma - Mat , = 1.2 6 is the a x- 
ial charge of the nucleon, = y/W^+w^ is 

the pion energy and UNNik), WAfA(fc), • • ■ are the 
NNtt, NAtt, . . . form factors associated with the 
emission of a pion of three- momentum k. The 



2 



form factors reflect the finite size of the baryonic 
source of the pion field and suppress the emission 
probabiUty at high virtual pion momentum. As a 
result, the self-energy integrals are not divergent. 

The leading non-analytic (LNA) contribution 
of these self-energy diagrams is associated with 
the infrared behavior of the corresponding inte- 
grals; i.e. the behavior as fc 0. As a conse- 
quence, the leading non-analytic behavior does 
not depend on the details of the form factors. 
Indeed, the well known results of chiral pertur- 
bation theory [^Q are reproduced even when the 
form factors are approximated by u{k) = 9{A—k). 

Of course, our concern with respect to lattice 
QCD is not so much the behavior as rriTr — > 0, but 
the extrapolation from high pion masses to the 
physical pion mass. In this context the branch 
point at = AM^, associated with transitions 
of — > A or A ^ iV, is at least as important as 
the LNA behavior near TOtt — 0. 

Heavy quark effective theory suggests that as 
— > c» the quarks become static and hadron 
masses become proportional to the quark mass. 
In this spirit, corrections are expected to be of or- 
der l/niq where ruq is the heavy quark mass. The 
presence of a cut-off associated with the form fac- 
tor acts to suppress the pion induced self energy 
for increasing pion masses, as evidenced by the 
in the denominators of Eqs. (||) , (^) and (||) . 
While some dependence in the form factor is 
expected, this is a second-order effect and does 
not alter the qualitative feature of the self-energy 
corrections tending to zero as 1 /m^ in the heavy 
quark limit. 

Rather than simplifying our expressions to just 
the LNA terms, we retain the complete expres- 
sions as they contain important physics that 
would be lost by making a simplification. We 
note that keeping the entire form is not in con- 
tradiction with xPT. However, as one proceeds to 
larger quark masses, differences between the full 
forms and the expressions in the chiral limit will 
become apparent, highlighting the importance of 
the branch point and the form factor reflecting 
the finite size of baryons. 

As a result of these considerations, we propose 
to use the analytic expressions for the self-energy 
integrals corresponding to a sharp cut-off in or- 



der to incorporate the correct LNA structure in 
a simple three-parameter description of the m,r 
dependence of the N and A masses. In the heavy 
quark limit hadron masses become proportional 
to the quark mass. Hence we can simulate a lin- 
ear dependence of the baryon masses on the quark 
mass, TTiq, in this region, by adding a term involv- 
ing m^. The functional form for the mass of the 
nucleon suggested by this analysis is then: 

Mn ^ aN + f^Nml + aNNi^N) + CTNAi^N) , (5) 

while that for the A is: 

Ma = OA + Pa^I + (taa(Aa) + ctawCAa) ■ (6) 

1.2. Model Dependence 

The use of a sharp cut-off, u{k) = 9{A — k), 
as a form factor may seem somewhat unfortu- 
nate given that phenomenology suggests a dipole 
form factor better approximates the axial-vector 
form factor. However, the sensitivity to such 
model-dependent issues is shown to be negligi- 
ble in Fig. ^. There, the self-energy contribu- 
tion aNN{= ctaa) for a 1 GeV dipole form fac- 
tor (solid curve) is compared with a sharp cut-off 
form factor combined with the standard a + jS 
terms of or (^. Optimizing A, a and /3 pro- 
vides the fine-dash curve of Fig. ^. Differences are 
at the few MeV level indicating negligible sensi- 
tivity to the actual analytic structure of the form 
factor. 

Here we have focused on the pion self-energy 
contribution to the N and A form factors. Only 
the pion displays a rapid mass dependence as the 
chiral limit is approached. Other mesons partici- 
pating in similar diagrams do not give rise to such 
rapidly changing behavior and can be accommo- 
dated in the a + (3 rn^ terms of (||) or (^) . More- 
over, the form factor suppresses the contributions 
from more massive intermediate states including 
multiple pion dressings. Other multi-loop pion 
contributions renormalize the vertex and hence 
we use the renormalized coupling gA as a mea- 
sure of the pion-nucleon coupling. 

2. ANALYSIS 

We consider two independent dynamical- 
fermion lattice simulations of the N and A 
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Figure 2. The self-energy contribution CTArAr for 
a 1 GeV dipole form factor (solid curve) is com- 
pared with a sharp cut-off form factor 0(Ajv — fc) 
(fine-dash curve). Self-energy contributions ctata 
(dot-dash) and ctaw (long-dash) for a 1 GcV 
dipole are also illustrated. 



masses. We select resuhs from CP-PACS's [| 
12^ X 32 and 16^ x 32 simulations at /3 = 1.9, and 
UKQCD's § 123 X 24 simulations at /3 = 5.2. 

Figure ^ displays fits of (||) to the lattice data. 
In order to perform fits in which A is uncon- 
strained, it is essential to have lattice simula- 
tions at light quark masses approaching ^0.1 
GeV2. 

It is common to see the use of the following 
xPT-motivated expression for the mass depen- 
dence of hadron masses, 



(7) 



The result of such a fit for the N is shown as 
the dashed curve in Fig. || The coefficient of 
the term in a three parameter fit is —0.761. 
This disagrees with the coefficient of —5.60 known 
from xPT (which is correctly incorporated in (|^) 
and illustrated as the solid and dash-dot curves 
of Fig. ^) by almost an order of magnitude. This 
clearly indicates the failings of (0). 

The dotted curve of Fig. || indicates the leading 
non-analytic term of the chiral expansion domi- 
nates from the chiral limit up to the branch point 
at = AM ~ 300 MeV, beyond which xPT 
breaks down. The curvature around m-^ — AM, 
neglected in previous extrapolations of the lattice 
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Figure 3. A comparison of phcnomcnological fit- 
ting functions for the mass of the nucleon. The 
solid curve corresponds to our preferred fit of the 
functional form of (||) with A constrained to re- 
produce a 1 GeV dipole form factor. The dash- 
dot curve illustrates the unconstrained fit. The 
three parameter fit (dashed) corresponds to let- 
ting 7 of (^) vary as an unconstrained fit param- 
eter. The dotted curve corresponds to (^ with 7 
set equal to the value known from xPT. The lat- 
tice data from are CP-PACS (solid) and UKQCD 
(open), each with a 5% scale change to provide 
consistency. 



data, leads to shifts in the extrapolated masses of 
the same order as the departure of lattice esti- 
mates from experimental measurements. 
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